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' I. INTRODUCTION 

, Quantum electrodynamics is known to be the best studied example of quantum field theory. Mainly, this is due 
£>• ■ to the weakness of the fine structure (coupling) constant, a w 1/137, what allows to perform many perturbative 
calculations as power series in a with an incredibly high accuracy. Despite the smallness of a, even in the realm of 
quantum electrodynamics, there are some questions that theory has not answered yet. In this paper, in particular, 
we address the problem of derivation of the low-energy effective action which at present is solved only partially for 
QED. 

' The low-energy effective action in quantum electrodynamics describes the dynamics of the electromagnetic field, 
Q\ I assuming that the production of the on shell fcrmions is absent or negligible. Apparently, such a description is self- 
consistent only if the fermions are massive and the characteristic photon energies are sufficiently small. The mentioned 
two conditions, as is clear, are necessary to suppress the process of the particle-antiparticle pair creation (on-shell). 
Intuitively, the low-energy effective theory is obtained from quantum electrodymanics by "integrating out" the 
' fermion field. After doing so, one arrives at a nonlinear theory that involves only the electromagnetic field degrees of 
freedom. In terms of the S-matrix language, one considers just those processes in QED which contain only photons 
among the asymptotic scattering states. The fermions, on the other hand, appear only through the internal loops by 
producing all kinds of photon vertices. 

The problem of deriving the effective action is an old one. Its roots go back to the well known papers of Heisenberg 
and Euler and Weisskopf [^j. There, for the first time, the effective action in QED (for the case of a constant 
electromagnetic field) was derived. From the viewpoint of application, the derived effective action contains, for 
example, the information on the photon-photon scattering at the tree level. It was this scattering process, in fact, 
that motivated consideration of the problem in Ref. [0J§], in the first place. Later, some further progress was achieved 
by Schwinger Q who, by using the proper time technique, rederived the result of Refs. [jl]j2| and, in addition, gave a 
nice interpretation to the imaginary part of the effective action in the case of a constant electric field. 

Obviously, the next most natural step in deriving the low-energy effective action in QED would be to take into 
account the effect of small deviations from the constant configuration of the field. In other words, the problem is to 
obtain the effective action as an expansion in powers of derivatives of the field strength. It turns out, however, that 
the latter is very difficult to accomplish (see for some early attempts in this direction) unlesss the weak field 
approximation is used. In this connection it is appropriate to mention that, in the weak field limit, the expansion is 
known up to four derivatives with respect to the field strength [pi. Our approach, on the other hand, does not involve 
any assumptions about the weakness of the background field. 
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A real progress in solving the problem started with the result of Ref. H] where an elaborated method, which, in 
principle, leads to a general result for the derivative expansion in QED, was presented. Because of the complicated 
character of the method, however, the explicit expression applicable to the most general case of the electromagnetic 
field background was not presented there. Recently, the derivative expansion of the effective action was obtained in 
the case of (2 + l)-dimensional QED This latter is a quite general result, containing all the terms quadratic in 
derivatives of the field strength with respect to the space-time coordinates. Finally, in our previous paper ||, we 
obtained a similar result for the effective action but in (3 + l)-dimensional QED. As in (2 + l)-dimensional case, it 
was given in a covariant form valid for the most general constant component of the electromagnetic field background 
what, as we will see later, is much more complicated problem than that in 2+1 dimensions. 

For completeness, we mention that some related interesting results were obtained in Ref. jl(j for QED and in 
Refs. [|ll]jl§ for non-Abelian gauge theories. 

In this paper we extend our method, which was originally presented for the case of (3 + l)-dimensional QED 
to QED in 2 + 1 dimensions. In particular, we obtain the derivative expansion of the effective action which includes 
up to two space-time derivatives of the electromagnetic field and, further, we formulate the Feynman rules for the 
perturbative expansion of the one-loop effective action in the number of derivatives. We also derive the explicit 
expressions for the derivative corrections to the imaginary part of the effective action in an external electric field. And 
finally, as a byproduct, we resolve the controversy posed in |]l3f where a result different from that of [^| was presented. 

The paper is organized as follows. In Sec. || we outline the general method developed in our previous paper 
Sec. |lj is devoted to solving some technical problems in dealing with functions of the matrix argument F^ v . Then, 
in Sees. |lV]and VII, we present the main res ults o f our paper, namely, the derivative expansions for spinor and scalar 
QED, respectively. In Sees. [y|-|v| and Sec. |VIIl[ p we calculate the derivative expansions for two particular cases 
of the external electromagnetic field, the purely magnetic and purely electric backgrounds, in both 2+1 and 3 + 1 
dimensions. Finally, in Sec. [x|, we develop the Feynman diagram technique for generating the perturbative expansion 
in the number of derivatives. Four appendices contain different formulas used throughout the main text. 



II. DERIVATIVE EXPANSION OF THE ONE-LOOP EFFECTIVE ACTION IN QED 

Let us start from the general formalism which was originally developed in [^) for (3 + l)-dimensional quantum 
electrodynamics. While doing so, we will notice that, to a great extent, the method does not depend on the dimension 
of the space-time. We will pay special attention to all those places where it does depend. 

In this paper we restrict ourselves to the one-loop effective action. This is the same approximation which was used 
by Schwinger || in the case of a constant external electromagnetic field. 

As is known, the one-loop effective action in QED reduces to computing the fermion determinant 

W {1) {A) = J d n xC {1) = -i\jxT>et{iV - to) = - % - InDet (l>£ + ^a^F^ + m' 

= - l -J d n x( X \tr\n(vl + ^F^+m 2 )\x). (1) 

Here V = 7^2?^ and the covariant derivative is = <9 M + ieA^. By definition, = i[ , y fl ,j l ,]/2 and tr refers to 
the spinor indices of the Dirac matrices 7^. States are the eigenstates of a self-conjugate coordinate operator x^. 
Throughout the paper we use the Minkowski metric, i.e., r]^ = (1,-1,-1) or rj^ v = (1, — 1, — 1, — 1), depending on 
the actual space-time dimension. And in both 2+1 and 3+1 dimensions, we work with the 4x4 representation of 
the Dirac 7-matrices. 

For calculating the effective action in Eq. (|l]), we employ a version of the so-called worldline (or string-inspired) 
formalism developed in |^-^6|. Such an approach to an ordinary field theory, based on the path integral over one- 
dimensional world lines, was extended to the evaluation of Feynman diagrams for Green functions in higher loop 
orders |l7|-|i{^|. It has demonstrated its power reproducing known theoretical results in QED while allowing one to 
invoke new technique to study the theory's behavior in strong couplingregime pC)[ ]. For some recent applications of 
the worldline formalism as well as for an extensive list of references see |2l],|2^]. Note, however, that our method differs 
from the one commonly used in the literature by a choice of the worldline propagators, and is closer in spirit to the 



method used in 17 2 



With use of the formal identity ln(if + to 2 ) = — L 00 exp[— %t(H + m 2 )]e?T/r for introducing the proper-time 
coordinate r, the effective Lagrangian can be represented through the diagonal matrix elements of the operator 
U(t) = exp(— irH), 
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C (1 \A) = ^ J ^L e - im2T tr(x\exp(-iTH)\x), (2) 
o 

where the second order differential operator H is given by 

H = -n„IF + ^F^(x), n M = —iDp. (3) 

The matrix elements (x| exp(— irH)\x) entering the right hand side of Eq. (||) may be interpreted as the matrix 
elements of the evolution operator of a spinning particle with r and H being the proper time and the Hamiltonian 
of the particle. The corresponding canonical momenta are P^'s which obey the commutation relations [x^, P u ] = 
and are defined by (x\P^\y) — —id^8{x — y) in coordinate representation. Following the standard approach [ B4j , we 
represent the transition amplitude {z\U[r)\y) between points a;(0) = y and x(t) = z in terms of a path integral over 
the real and Grassmann coordinates, x^{t) and ip^it), as 

tr(z\U(r)\y) = N- 1 J V[x(t), ^(t)} exp jz J dt [L bos (x(t)) + L fer (t/;(t), x{t))] 1 , (4) 

where TV is a normalization factor, and 

1 dx dx u dx v 



i dip u 

Lf er (ip,x) = -ji>v-T t ie^ u ip F vX {x). 



(6) 



The integration in Eq. (Q) goes over trajectories x M (t) and ip^{t) parameterized by t <E [0, r]. The definition of the 
integration measure assumes the following boundary conditions 



x(0) = y, x(t) — z, "0(0) = —iP(t). 
We choose a special gauge condition for the vector potential A^(x), namely the Fock-Schwinger gauge [^5| 



{x v -y v )A v {x) = 0, 



(7) 



(8) 



which leads to the series 



M*) = \i* x - v x )Fxu{y) + \(x x - y x )(x° - y°)d„F x „{y) 

+ ±(x x - y x )(x° - y°){x» - y^d^Fx^y) + ... 
(x x - y x )(x^ - y" 1 ) . . . (x"« - y^) 

n=0 



E 



n\(n + 2) 



■ du n F\ v (y). 



(9) 



This choice of the gauge for the vector potential turns out to be very convenient for developing a perturbative theory 
in the number of the derivatives of the electromagnetic field with respect to the space-time coordinates. 

Carrying out the change of the variable x(t) for x'(t) = x(t) — y in the path integral in Eq. (Q) (henceforth we omit 
the prime) and substituting Eq. (0) into Eq. ([I]), we obtain 



tr{z\U(r)\y) — N^ 1 / D[x(t),^(t)]exp 



dtl 



iff -r XF ^ 



x exp 



(10) 



The new boundary conditions for x(t) are a;(0) = and x(r) = z — y. Notice, that F^ u in Eq. ( ]To| ) does not depend 
on x(t). As follows from Eqs. (||), (^]) and (g), the expressions for the interacting terms, L l f^[x) and L l ^ r (x,%l)), 
containing derivatives of F^ with respect to coordinates, take the form 
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j int 
^bos 



n!(n + 2) 



-x Ul (*)... x 1 ^ 1 (t) 



n=l 

e dx" A CT e dx v x a K 

Ti^uX^a 77~ X X "i~ uX.KTK 77~ X X X -\- 

3 at 8 at 



Lf: r { X ^) = -Y^-^Fx^ l ... v ^\t)r{t)x^{t)...x^(t) 



Here we use the conventional notation for the partial derivatives 



(11) 



(12) 



(13) 



Now we see that the problem of obtaining the derivative expansion reduces to the evaluation of the path integral in 
Eq. (O) in the framework of the perturbative theory with an infinite number of interacting terms given in Eqs. ( p"l| ) 
and (p_2[). Fortunately, for computing the effective action that includes only a finite number of the derivatives, it is 
sufficient to consider only a finite number of the interacting terms. Later, we shall restrict ourselves to obtaining only 
the two-derivative terms in the action. So far, we continue developing the scheme for the most general case. 

As usual, introducing real and Grassmann external sources, the matrix elements of the evolution operator can be 
represented as follows 



tr(z\U(T)\y) = cxp {i / dt 



LI 



1 S 



i Sr)(t) 



L 



fer 



1 6 



i6 V (ty sat) 



x Z T [ri,i]{z;y) 

where the generating functional is just the Gaussian path integral, 



(14) 



Z T [ V ^}(z:y) = N- 1 J D[x{t) , ^(t)] exp 



1 dx v dx u 



dx v 



2 dt dt 



dt 



x exp 



(15) 



The calculation of this generating functional reduces to obtaining the "classical" trajectories for x v (t) and ?/v(i), 
satisfying the appropriate boundary conditions, and to computing the determinants of the one-dimensional differential 
operators, 



and 



2 = irin V — - 2ieF fiU , 



(16) 



defined on the interval [0, r] with the periodic and antiperiodic boundary conditions for their eigenstates, respectively. 

The "classical" trajectories are easily obtained by solving the equations of motion that the bosonic and Grassmanian 
worldline actions in Eq. ( |l5| ) require. So, we arrive at 



„2eFt 



(z - y) v + 



dt' 



juFt _ ! (e 2eF ( T -*') - l) f (e 2 ^( t -*') - l) 



eF 



9(t-t')- 



eF 



(17) 



and 



r cl (t) 



^ (^ F{t - t,] (o(t-t')- TT ±^yut' ) . 



(18) 
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Then, the result of the path integration in Eq. (na) for the case of the coincident arguments z — y = x reads 



Det(0 2 ) ( i Qhos ,, l Q}er , 



Z T [ V , x) = C o y exp l^-STM ~ ^tilt] ) > (19) 

where the normalization constant Co should be determined by comparing the result with the Schwinger's one, or by 
satisfying the normalization condition 

Z T = [r),£](z;y) = 6(z-y), (20) 

which is equivalent to the operator equality U(0) = 1. The prime in Eq. ( |l^ ) denotes skipping a zero mode in the 
definition of the determinant. With our normalization convention for the determinants (see the next section), it is 
easy to check that the overall factor Co = — i/(27rr) 2 in 3 + 1 dimensions and Co = exp[— z7r/4]/[2(7rr) 3 / 2 ] in 2 + 1 
dimensions. 

The expressions for S^? a and S[f r are quadratic forms in the external sources 

T T 

S b c t s M = j 'dtxj dt 2 vAh)D v x (t 1 ,h)r, x (t 2 ), (21) 



T T 

sir® = j 'dtxj dhuti)s v x (tuh)e(t 2 ), (22) 



where the Green functions are given in terms of functions of the matrix argument F^ v 



x ' ; 2eF 



e(h - t 2 ) (l - e 2eF ^~^) + coth(eFr) (l + e 2eF ^'^) 



e eF(r~2t 2 ) + e eF(2ti-r) 



sinh(eFr) 



(23) 



S(t u t 2 ) = - [ e (t! - t 2 ) - tanh(eFr)] e 2eF ^\ (24) 

Substitution of Eqs. @, @ and @ into Eq. @) leads to the expression for tr(x\U\x). After expanding the 
exponent in powers of the operator valued interacting terms, L^ s and L l p^ r (containing functional derivatives with 
respect to the sources rj^t) and ^(t)), one has to calculate the result of the derivative action on the generating 
functional. Starting from this point, we have to restrict ourselves to a specific finite number of the derivatives in the 
effective action. As we mentioned before, in this paper we are interested in the two-derivative terms (see Sec. [x] for 
some discussions on computing the higher order approximations). Therefore, we obtain 

/ T T T 

tr(x\U{ T )\x)= (l + i f dt[V 2 (t) + W 2 (t)]-^ [ fdhdt 2 [^(fOViCtaJ + WkfajWiCfe)] 



x 

dhdhV^Wxih) J Z r [ri,£](x, x) 







(25) 

77=04=0 



where, as follows from Eqs. (|l l|) , ( |12| ) and (|14|), the vertex generating operators are 

, x i d S 3 

Vi(t) = —eF v \ t n lira 



3 u ^t a ^t dt ^(t )%(t)%(*) 

1 d S 4 

V 2 [t) = -eF v \ tliK lira 



W^t) = —e.F v \ ttl 



8 ' ^ to-t dt 5T} v {t )5r]x{t)5r)^{t)5r] K {t) ' 
5 2 S 



i 



sum*® 

s 2 s 2 



W 2 (t) = -eF^ K (26) 
2 i^W^AW or}^{t)6r} K {t) 
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Substituting the generating functional ( |19| ) which depends on the Green functions ( |23| ) and (|24]) , we rewrite Eq. (|2£ 
in the form 



*r<a;|I7(T)las> = C 1 



Det(0 2 ) 



1 - oeF„ A „ 



dt 



Det'(Oi) 

u 

+D vp (t, t)D XK (t, t) + D VK (t, t)D x »(t, t) + 4S" x (t, t)D» K (t, t) 
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dtidt-2 



9£> w (l, 2) (S KCT (2, 2)S A,y (l, 1) - 2S kA (2, 1)S™(2, 1)) 

+6S"" t (2, 2) (z)" A (l, 1)£> W (1, 2) + 1)13^(1, 2) + 2)L> A ^(1, 1)) 

+ZT A (1, 1)£) CTK (2, 2)D W (1, 2) + 2D uX {l, 1) (l)^, 2)Z^ K (1, 2) + D a »(2, l)D Kp {2, 2) 
+1)^(1, l)D ap (2, 2)D Xk {\, 2) + 2IT K (1, 2) (Z) CTp (2, 2)L> A ^(1, 1) + D ap (2, l)L> Ap (l, 2) 
+ZT K (1, 2)-D ctA (2, l)D pp {\, 2) + D vp (l, 2)D ap (2, 1)D Xk (1, 2) 

+D" a (l, 2) (D xp {l, l)D Kp (2, 2) + D Xk {1, 2)D pp (l, 2) + D Xp (l, 2)D pfi {\, 2)) 
Here the dotted functions are defined by the expressions 

ir(i,2)Wlir(ti,i J ) 1 

Oti 

^(M)^lim|^(f„,t). 

to-»t OTo 



(27) 

(28) 
(29) 
(30) 



Having the representation (27) together with the Green functions fl23|) and (|24|), one is left with a need to perform 
the integrations over the proper time. This latter, however, may look like a rather complicated problem due to the 
necessity to disentangle the Lorentz indices while doing the integration. In the next section, we show how this problem 
can be solved. 



III. HOW TO DEAL WITH FUNCTIONS OF MATRIX ARGUMENT 

In the previous section we developed the general method for calculation the derivative expansion in QED. However, 
there was not given an explicit final expression, since we needed a technique dealing with functions of the matrix 
argument F^. Below we show, following the method of ]2q] , how to deal with those functions as well as how to 
calculate the determinants of the differential operators in Eq. (|l6|) . 

Let us begin by introducing notations that we are going to use below. When working with the electromagnetic 
field strength tensor, it is usually very convenient to introduce the invariants built of the field strength. In (3 + 1)- 
dimensional theory, the standard choice of the two independent invariants reads 

T = --F^F^, G = \ e pvXK F Xfi F^. (31) 
4 8 

In our calculations, though, it will be more convenient to work with the following couple of invariants 

K+ = yj V^ 2 + Q 2 K- = yj V^ 2 + G 2 (32) 

As for the (2 + l)-dimensional theory, there exists only one independent invariant built of the electromagnetic field 
strength, and it is given by the expression analogues to T in Eq. (|3l|). 

Now we proceed to the case of (3 + l)-dimensional QED. It is this case that was considered in p6fl . The authors of 
that paper introduced the set of matrices A v , X s with j 6 {1,2,3,4}, 



G 



.4 



2(/ 7 2 -/7) : 



(33) 



where 



f 1<2 =±iK_, h. A = ±K+- 
fi. 2 = TK + , f 3 . 4 = TiK-. 



(34) 
(35) 

The main property of the matrices ([33]) that we are interested in are their (left and right) contractions with the 
field strength tensor, 



F v A {i)X ^ - A^F Kfi = fiA^^ 
Other useful properties of these matrices that will be used below are 



fJLV 



(36) 



(37) 



As follows from the property in Eq. (36), for any function <&(F) of the tensor argument F^, we get 



(38) 



Matrices with similar properties can also be introduced for (2 + l)-dimensional tensor F^ as well. Indeed, the 
following set of matrices 



a^ v - _ 

(±D ~ 2 V 2T 



1 f {F 2 ) tJ - v F^" \ 



IT 



(39) 



in the (2 + l)-dimensional case have properties similar to those in Eqs. (36) and ([37|). As is easy to check directly, 
their eigenvalues are 



f ±1 = ±V2F, f a = 0. 



(40) 

In particular, for the Green functions (^) and (|24|), which are functions of the tensor argument F^, we obtain the 
following representations, 



2efi 

e e/ 3 (r-2t 2 ) + e e/ J (2t 1 -T) 

smh(efjT) 



_ e 2e/ J -(t 1 -t 2 ) 



■coth(e/j-T) (l + e 2 ^^ 1 -* 2 )) 



(41) 



S" A (ti,ta) =J2 A U)l - t2 ) - tanh(e/ i r))exp[2e/ 7 (t 1 - t 2 )]. (42) 

j 

As is seen, in the case of vanishing field, the propagators D 1/A (t 1 , i 2 ) and S vX (ti,t 2 ) coincide with those used in Jl7| , p3]| . 

Another problem is related to calculating the determinants of the operators ([!(]). The latter are nothing else but 
products of all eigenvalues of the operators. Once again, making use of the matrices in Eq. ( ]33| ) or in Eq. (39) for 
(3 + 1)- or (2 + l)-dimensional cases, respectively, we look for the eigenvectors of the operators 0\ and 2 in the form 

4)W=^(i)A«VW (43) 
V U) (t)=Al j)x er,{tl (44) 

where a A and £ A are constant nonzero vectors, cf> and r\ are scalar functions of t. As a result, the problem of 
obtaining eigenvalues reduces to solving ordinary differential equations for the scalar functions <p and rj with appropriate 
boundary conditions. 
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Now, it is easy to check that, up to an unimportant constant, the corresponding determinants read (note that we 
skip a zero mode of the operator Oi) 



De^ 3+1 )(0 1 ) 



sinh 2 (erK + ) sin 2 (erK_) 



{erK + ) 2 {erK-) 2 ' 
L>ei (3+1) (0 2 ) = cosh 2 (erK + ) cos 2 (erKJ) 



in the case of QED in 3 + 1 dimensions, and 



Det^ 2+1 \0 1 ) = 



sinh 2 (er\/ir) 

2l 



(45) 
(46) 

(47) 
(48) 



Det {2+1) {0 2 ) = cosh^(erV2J : ') 
in the case of QED in 2 + 1 dimensions. To obtain these results we used the following formulas for infinite products 



n i 



sinh : 



n i 



Ax 2 



7T 2 (2n+ 1)' 



= cosh ; 



(49) 



and similar ones with replacement x — > iy. 



IV. GENERAL RESULT IN SPINOR QED 

By making use of the results from the previous section, we can proceed with the calculation of (p7|). 

After substituting the Green functions ( pd| ) and (f42|), as well as the explicit expressions for the determinants of the 
operators 0\ and O2, a straightforward, though tedious computation gives the result for the diagonal matrix element 
of the U(t), 

tr{x\U{r)\x) =tr(x\U(T)\x) 



-^e 2 F^,F aK , p J] 9CF w (fj,fi,fk)A l fiA%A% ) + 9C^ w (f j} f h f k )A$A%A% 
j,i,k \ 

+6CY w (f j ,f l J k )A^ ) (^A^j +^^,) +60^,^/,,)^^ 
-CY V (fj, fu fk) (AftAtfAft + A%A%Afc + 2A^A^ h 
-C% v UiJu fk) (A^AftA^ + AftAftAfc + 2A»$A«*A* X 



-2G 



-a 



vv 



fufk) (A$A%A% + A%A$A$ - CX v (f J , fuh)A^A 

(f. f, f,\A UK ( A Xg A^ p A- A Xp 4 A " T 
Uj, JhJk)A(j) [ A (i) A (k) +A (i) A (k) 



op 

(k) 



(50) 



where the explicit expressions for the coefficients C* Y (a, (3, 7) (with X, Y £ {V, W}) are given in Appendix [A| and 
the diagonal matrix elements tr{x\U(r)\x)o correspond to the non-derivative case, 



tr{x\U{r)\x)^ +1) =- 



An 2 i 



:(eTK_)(eTK + ) cot(erA'_) coih(eTK + ) 



in 3 + 1 dimensions, and 



tr(x\U (r)\x)^ = ^J^ ] (erVW) coth( e rV2^) 



(51) 



(52) 
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in 2 + 1 dimensions. 

Equation ( [sof ) (along with a similar one for scalar QED) is the main result of our paper. Note that the renor- 
malization of the effective action (||) formally reduces to (i) performing a subtraction (precisely the same as in the 
original Schwinger's paper ra) of a term containing no derivatives of field strength with respect to coordinates, and 
(ii) changing all bare quantities for the renormalized ones, e — * en and — ► A^, defined as follows: 



z 3 e, 



A, 



,-1/2 



A, 



1 + Ce 2 



(53) 



where 



oo 



q(2 + 1) = 



1/A 2 
oo 

i r ds 







exp(— sm 2 ) 



Qirm 



(54) 
(55) 



and A is an ultraviolet cutoff in (3 + l)-dimcnsional QED. 

After subtraction and conversion to the renormalized quantities the effective action becomes finite in the limit 
A — > oo. Since the derivative part of the effective action depends on e and A^ only through the product eA^ = e^A^ 
it does not change its form and no further renormalization is required to make the derivative part well defined (below 
we use only renormalized quantities, although we always omit the script "R" in their notation). 

By using the asymptotic behavior of the coefficient functions (given in Appendix |a|), one easily finds the following 
expansion of tr(x\U(r)\x) in powers of t, 



tr{x\U{f)\x 

1 



tr(x\U{T)\x) 



20 



180 



2 



p vX - p n 



(56) 



As is clear, this is the weak field limit of our general result in spinor QED. In the effective action, the given order in 
t results in the two-derivative corrections of the order 1/rn 2 : 



V™ 2 ~ 7207TTO 2 

in 3 + 1 dimensions, and of the order 1/m 3 



fpvX,np 



(57) 



C 



(2+l)spm 
1/m 3 



720m 3 



9P" A ' p M 



(58) 



in 2 + 1 dimensions. 

The expansion in Eq. (|56|) was obtained earlier in the heat kernel approach . While the latter is a perfect tool for 
deriving the effective action in the weak field limit, it is not very useful when the field becomes strong. Our approach 
here, on the other hand, is free from such a limitation and the general result in Eq. ( |50| ) contains all the two derivative 
terms like dFdF(F/m 2 ) n where n is an arbitrary positive integer and the Lorentz indices (not shown) are contracted 
in al l po ssible ways. To substantiate this claim, we present the next to leading terms of the weak field expansion in 
Eq. Q in Appendix |b| 

As we saw above, the formal expansion in r corresponds to an expansion of the effective action in the inverse powers 
of the mass parameter. This means that, while making use of such an expansion, one cannot get any reliable results 
in the limit of the vanishing fermion mass. This, in particular, is the main reason why the authors of [fl3|| , who used 
an expression like (|5^), came to a wrong conclusion about the absence of corrections to the one-loop effective action 
coming from inhomogeneities of a static magnetic field when m — * 0. Such a conclusion "contradicts" the result of 
Ref. pi. The latter, as we will see, completely agrees with our result for the derivative expansion. 



V. SPINOR QED IN 2+1 DIMENSIONS 

Let us consider the case of the purely magnetic field background to which a special attention was paid in ||. To 
proceed with analyzing this case, note that the electromagnetic field strength tensor takes the following form, 
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F^(x) = B(x)F^, 



(59) 



where B(x) is a pseudoscalar function coinciding with the magnetic field strength and is a constant matrix 
with the only nonzero components F 12 = — F 21 = 1. As is seen it satisfies the following normalization condition: 



F""F MI/ = 2. 



To reduce the general result presented in Eq (p0|) for the particular choice of the field given in Eq. fl59|), we have 
to use the properties of A^'s presented in Sec. ID. Just to get feeling how they work, let us consider an example, 



3,1 



d ^^C w (f J J l )f u) A^ ) 



3d 



B 



C 



w 



(60) 



In this derivation, we made use of the Bianchi identity. We recall that the latter should be satisfied since the 
electromagnetic field was introduced in the theory through the vector potential by minimal coupling. The identity 
itself reads A^d„B = 0. The direct consequence of it is the independence of the magnetic field, for the particular 
choice (p9|), on the time coordinate. By noticing that the matrix A^y as well as any other from the set, does not 
depend on B(x) we obtain the secondary identity, A^d^d v B = 0, by differentiating the original one. It is this last 
form of the Bianchi identity that was actually used in our derivation in Eq. (|60| ) . 

The other expressions, similar to that in Eq. (|60|), along with the functions like C w (\ / 2J 7 , \J2T~) are listed in 
Appendix ^|. 

The final result for the derivative part of the diagonal matrix element (|50|), for the particular choice of the field 
configuration in Eq. ([59]), reads 



tr{x\U{r)\xY d 



(2+1) 



^e 2 {d,B) 



1 



(47r|eB|) 3 / 2 



to 



(3w 2 F 4 - 3wF 3 - Auj 2 Y 2 + 3ojY 



- OJ 



ie 2 {diBY 



(47r|eB|) 3 / 2 2 du; 3 



(u) coth w) . 



(61) 



where to = ir\eB\, Y — cothw, and (diB) = diBdiB. Substituting the last expression into Eq. (|2|), we come to the 

integral representation for the derivative part of the effective Lagrangian (we perform the change of the integration 
variable r for w = ir\eB\), 



->{2+l)spin 



(B) 



(diB) 2 



4(47r|eB|) 3 / 2 J ^ 



dui 



■. exp 



eB 



dco 3 



(u> cothw) . 



(62) 



The last expression coincides with the result presented in g (note that in notation of || diBdiB = AdBdB). One can 
be convinced that the integrand in (|6^) is a negative function what means that inhomogeneities of the magnetic field 
background, in approximation under consideration (one- loop and two derivatives), lead to the reduction of vacuum 
energy density for any value of the ratio m 2 /\eB\. The latter situation does not, however, prove that a spontaneous 
generation of a non-homogeneous magnetic field happens in QED since the sign of the two derivative term in the 
expansion of the effective action is not a sufficient argument for making a conclusion of that kind p8[ . 

We would like also to give another representation for the derivative part of t he L agrangian in terms of special 
functions. To get it, we need to perform the integration in (62) by parts (see Eq. (D5) in Appendix ^). Here is such 
a representation, 



C 



(2+l)spm 
der 



(B) = - 



(diB) 2 



\/27r(4|e£|) 3 / 2 



2\eB\ 



5C 



C( 5 ,i 



3 m 2 
+ 2kB| 



2 eB 



ci-ii 



2\eB\ 



2\eB\ 



2\eB\ 



2eB 



(63) 
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Often, in the limit of large or small values of the external field, it is more convenient to work with the asymptotic 
expansions of the effective action rather than the exact expression as in Eq. ([33]). First, let us consider the case 
to 2 <C \eB\. Then, using the last representation, we easily derive the following asymptotic expansion, 



C 



(2+l)spm 
der 



(B) 



e 2 (d l Bf 



V2(4n\eB\)V* ^ 



5-2fc / 1\ / 3 
E I k + n ) C( k - 



k\ 



2\eB\ 



(64) 



In order to get the asymptotic expansion for m 2 3> \eB\, we make use of the integral representation in Eq. ( p2[ ) and 
obtain 



r {2+l) spin , 



e\d % B) 2 ™ B 



27r 3 / 2 m 3 ^ (2fe + l)! 



2fe+4 r(2fc+ 3\ f2\eB\ 



2/,- 



(65) 



where B^ are the Bernoulli numbers. 

Now, let us consider the case of the purely electric field background. Without loosing the generality, we assume 
that the field is directed along the first axis of the two-dimensional space. Again the field strength tensor is factored 
similar to (B9l), 



F fu/ (x) = E(x)F^ 



(66) 



where E(x) is the magnitude of the electric field. Now the constant matrix F^ v has nonzero components F 10 = — F 01 = 
1, and satisfies the normalization condition: F^F^j, = —2. The general expression ( [sol) simplifies considerably for 
our choice of the background field. And the derivative part of that expression now reads 



tr{x\U{r)\x) { ?l x \E) 



iexp(-i7r/4) e 2 (d\\E) 2 
(4n\eE\) 3 / 2 

«exp(-i7r/4) e 2 (d\\E 
~ {4i:\eE\) 3 / 2 2~ 



(3w 2 F 4 - iuoY 3 - Auj 2 Y 2 + iujY 
) 2 d 3 

~V^~r^ (wcothw) , 
aiv 3 



(67) 



where now u> = r|eJ5|, Y = cothw, and (d\\E^j = {d^Ed^E — d\EdiE) . Here we used the Bianchi identity again 
to show that the electric field does not depend on the second spatial coordinate. Substituting this expression into 
Eq. (0), we come to the integral representation for the derivative part of the effective Lagrangian, 



^»(2+l)spm 

der 



(E) 



exp(-i7r/4)e 2 (d\\E) 
4(47r|e£|) 3 / 2 



: exp 



eE 



d 3 

—-^r (u> cothw) . 
dio 3 



(68) 



As expected in the case of an electric field background, this derivative correction to the effective action contains a 
nonzero imaginary contribution. A convenient representation of the latter can be obtained in the following way. First, 
in Eq. (|6"8|), we switch to a new variable, z — iu>, so that the integration runs along the imaginary axis of z from zero to 
ioo. Then, we move the integration contour to the real axis of z. As is easy to check, the integrand has poles at z = nn 
(n =1,2,.. .). As a result, the real and the imaginary contributions get naturally separated. Indeed, the real part of 



,(2+l)spm 
J der 



is given by the principal value of the integral along the Re(z) axis, while the imaginary part appears due 
to the integration along the infinite set of the vanishingly small semi-circles above the poles, z = Tin + eexp[i(7r — <f>)] 
(where < <f> < ir and e — * at the end). In this way, we easily obtain the imaginary part of the right hand side in 

Eq. m, 



TmC 



(2+1) spin 



(E) 



1 



2 s TT 3 \eE\ 3 / 2 ^ n 5 / 2 

1 1 n— 1 



15 + 18- 



eE 



12 



exp 



eE 



irm 2 n 



2 \ 3 ' 

■Km n 



eE 



(69) 



We note that the result of the summation in the last expression (as well as in similar formulas later on) can be given 
in terms of the polylogarithmic function Li^ (x) . Equation ( |69|) determines the correction to the probability of the 
particle-antiparticle pair creation (by definition, the probability density is W = 2XmC) in an external electric field 
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due to small inhomogeneities in space-time. We emphasize that the correction due to a time derivative of the field 
has the "wrong" sign, i.e. it works against the particle creation. The gradient in the space direction parallel to the 
field strength, on the other hand, amplifies the process. 

As is known, in the case of constant electric field, the imaginary part of the effective Lagrangian is given by 

This as well as the first correction due to the derivatives remain finite even in the limit of zero fermion mass. Despite 
of this fact, we still expect that the derivative expansion (with the electric field background) may fail in the limit of 
vanishingly small mass due to higher orders in the number of derivatives. Below we shall see that the same is true in 
the spinor QED in 3 + 1 dimensions as well. 



VI. SPINOR QED IN 3+1 DIMENSIONS 



As was mentioned at the beginning of the paper, the derivative expansion in QED4 was also studied in |7j. The 
result of that paper was presented in an explicit form for the special class of the electromagnetic field configurations, 



0. 



F" u (x) = <f>(x)F» v , 



(71) 



where $(x) is a slowly varying function that defines the magnitude of the field, and is a constant matrix. For 
convenience, let us normalize the matrix F Ml/ by the condition: F^F^,, = 2. Then the scalar function <!>(x) is nothing 
else but y/(-2T). As was shown in our previous paper [B| , the general result for the diagonal matrix element (pG) in 
the case of field ( f7l| ) reduces to the same result as was presented in |7j , 



tr{x\U(r)\x)^ 1) ^) = 



(4tt) 2 t $ 2 



ZujY 3 



4lu 2 Y 2 + ZujY 



(47T) 



d u $d>*<f>w d 3 , 

(Ld COthw) , 



<J> 2 2 duJ 3 



(72) 



where u> = re<l>, Y = cothw. As in the (2 + l)-dimensional theory, here we used the Bianchi identity, which this time 
reads 



(if + (F 2 ) 



d u $ = 0. 



(73) 



In the case of magnetic field along the third axis, for example, this condition means that the specified field cannot 
depend on the time and the third spatial coordinates, while in the case of electric field along the first axis, it cannot 
depend on the second and third spatial coordinates. 

Now, let us consider two particular cases of external field that we studied in 2 + 1 dimensions: purely magnetic and 
purely electric field backgrounds. Both of them are just different possibilities of that given in Eq. ([FT). 

Thus, in the case of magnetic field (along the third axis in space) we come to the following integra representation 
for the derivative part of the effective Lagrangian, 



(3+l)spin 



(diB) 



B7r) 2 |eB| 



exp 



eB 



d 3 
dtu 3 



(u cothw) 



(74) 



Resembling the situation in 2 + 1 dimensions, inhomogeneities of the external magnetic field tend to reduce vacuum 
energy density for any value of the ratio m 2 /\eB\. 

Performing integration in the right hand side of Eq. ( [74] ) by parts (see Eq. (D6) in Appendix we find the 
following representation (for the representation of the part of the effective action without derivatives in terms of 
special functions, see fl30|), 



.(3+l)spin 
J der 



(B) 



e 2 (d. t B) 2 
(87r)2|eB| 

24C (-2,1 
G 



rn 



11 

~6~ \ \eB 

o 



eB 



1 m 



3 eB 



eB 



i> 1 + 



2eB 



\eB 



2\eB\ 
lnT I I 



- 24- 



eB 



C -1,1 



2\eB\ 



2\eB\ 



lnV27r 



(75) 
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As m 2 <C \eB\, this expression allows the following asymptotic expansion, 



r (3+1) spin, r,\ 



e 2 {d % Bf 



(8?r) 2 |eB 
m 



24C'(-2) + 



2m 1 



m 



2|eB|*^jfc + 4 



m ~ 



2\eB\ 



where C'(-2) ~ -0.030. As m 2 > \eB\, on the other hand, we obtain 



■*(3-\-l)spin 



(B) ~ — 



(2tt) 2 to 2 



2fc + 1 



fc=0 



2|eS| 



m- 



2 A- 



77? 



3|eE| 2|e5| 2 3|eS| 3 
2 x fc" 



(76) 



(77) 



In case of the electric field along the first axis, on the other hand, we obtain the following expression for the 
derivative part of the effective action, 



r (3+l)spin , px 
L der ( K ) 



ie' 



(d\\E) 



(87r) 2 |e£7| J to 



du> 

— exp 



eE 



d 3 

uj I — — r (u> cothw) . 
' auj 6 



(78) 



This expression has both real and imaginary part, as always happens in the case of an external electric field. Another 
representation for it is obtained by analytical continuation of ( f75| ) according to the rule \eB\ —> — i\eE\. The imaginary 
part though is easily extracted from (|78|) in a standard way, 



XraL 



(3+l)spir 
der 



(E) = 



(d«E 



oo 1 

T — 



2 6 n 4 \eE\ ^ 

' ' n— 1 



6 + 6- 



exp 



eE 



eE 



2 \ 2 

Km n 



eE 



2 x 3 ' 

■Km n 



eE 



(79) 



which determines a correction to the Schwinger result |^] for the imaginary part of the effective action in a constant 
electric field, 



Tmrt 3+ V spin (E) = Y -i exp 

Sir 6 n z 

n— 1 



eE 



(eE) 2 



exp 



eE 



(80) 



The result in Eq. ([79|) is in agreement with that of ■ 

As is easy to establish, both the Schwinger result for a constant field and the first correction due to derivatives are 
finite in the limit of the vanishing fermion mass. As we argued in the case of the (2 + l)-dimensional spinor QED, 
this may not be the case in higher orders of the perturbative expansion in the number of derivatives. 



VII. GENERAL RESULT IN SCALAR QED 

Now turning to the calculation of the derivative expansion for the scalar electrodynamics, one does not need to 
repeat all the calculations similar to those done in Sec. In order to see this, we recall that the effective one-loop 
Lagrangian in this case reads 

oo 

L {l)scal {x) = _. / *I( a .| I 7 6M ( T )| a .) e -im S r_ (81) 



The evolution connected with the transition amplitude, (z\Ub os (T)\y) , is described now by the Hamiltonian (compare 
with Eqs. (|) and ©) 

H bos = -II^IP, n„ = -ify + eA^x). (82) 
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Thus, omitting all terms originating from the fermion part in the expression (^), i.e. putting L l J^ r = in Eqs. (|Tc|), 
( |l4| ) and S[^ r = in Eq. ( |l9] ) , we come to the following expression 

(x\U bos (T)\x) = (x\U bos {T)\x) 



18 



-e z F vX .uF, 



v r vv (f- f, f, M 4 I/A 4 K<T 4 W 4- A ufl A Kp 4 A<T 4- 9 4"* 4 Kp 4 P '' T 

x 2^ °i \h'JhJk)[ A (j) A (i) A (k) + A (j) A (i) + ^(0 A (fe) 



+2Cr v (f 3 , fiJk) (AftAtfAft + A^A^j + CrUi,fu hWJ^A^ 
+cY v ifs, fuh)A v & + 



(83) 



The coefhcients used here are the same as in Eq. (pOj) . As for the non-derivative factors, they have the standard form, 

,(3+i) _ i (eTK_)(eTK + ) 



(x\U bos (T)\x) { 



in 3 + 1 dimensions, and 



(x\Ubos(T)\x)^ 



(Akt) 2 sm(erK ) sinh(erfC + ) 



(2+1) _ exp(-i7r/4) (erVZF) 



(4ttt) 3 / 2 sinhierVW) 



(84) 



(85) 



in 2 + 1 dimensions, as can be easily checked by using the expressions for the determinants given in Sec. Ill and by 
taking into account the fact that, because of spin degrees of freedom, we had the additional factor 4 for fermions. 

In the case of scalar theory, the renormalization of the electromagnetic field and charge is given by the same formulas 
( [53| ) but this time the corresponding constants read 



00 

C ' 3+1, = I^ / 7-<— '»■ 

1/A 2 

00 

C^+D = —L- [ *L cxp(- S m 2 ) = — . 



(86) 



(87) 



To get a result of the type as in B, one has to expand the coefficient functions in powers of proper time. Thus the 
expansion for (x\U bos (T)\x) (weak field limit) reads 



(x\U bos (T)\x) = (x\U bos {T)\x} 



,2 3 

1 - — F uA F 

30 vX,n 



180 



p" A - F ^ 



This expansion up to the order t 5 is given in Eq. ( |B26| ) in Appendix [b]. As in the spinor QED, it is useful only in 
the case of heavy scalar particles (weak fields), when the mass scale is much larger than all other scales in the theory. 

In the effective action of scalar QED, the expansion in Eq. ( |S8l ) corresponds to the following leading two derivative 
terms 



C 



(3+l)scaZ 
1/m 2 



in 3 + 1 dimensions, and 



C 



(2+l)scal 



1/r, 



720mn 2 



a 

720m 3 



RP^p ' 1 + 4P" a ,mp . + F yX F V 



(89) 



(90) 



in 2 + 1 dimensions. 
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VIII. SCALAR QED IN 2+1 DIMENSIONS 



Let us start by considering the case of an external magnetic field as in Eq. (p9|) ■ This time the derivative part of 
the general expression (R3) reduces to 



(x\U bos (r)\x)l^ 



ie 2 (diB) 



4(47r|eS|) 3 / 2 sinh( 



-(3ojY 3 - 3Y 2 - 2uY + 1) 



ie: 



(d t B) 2 



A— 

A(Air\eB\f/ 2 2 \doj z duo) Vsinh. 



(91) 



where uo — ir\eB\, Y — cothw. After substituting the last expression into Eq. (pjj), we come to the integral rep- 
resentation for the derivative part of the effective Lagrangian (after performing the change of integration variable 
r — > uo — ir\eB\), 



r (2+i) seal , e 2 (diB) 

L der ( B ) = 



(167r|eB|) 3 / 2 



: exp 



eB 



d 3 d_ 

duo 3 duo 



sinh uo 



(92) 



which coincides with the result presented in As in the case of spinor QED, there exists another representation of 
( |92| ) given in terms of special functions (see Eq. ( Dll ) in Appendix |5]) , 



,(2+l)scal 



(B) 



e 2 (diBY 



V2n(l6\eB\) 3 / 2 

2 x ^ 



20C 



3 1 



2' 2 2\eB\ 



- 18- 



eB 



c.-i i 



2' 2 2\eB\ 



eB 



c 



1 1 



2' 2 2 eB 



eB 



eB 



3 1 



2' 2 2eB 



(93) 



Numerical study of the integral in (|92|) shows that inhomogeneities of magnetic field background, in approximation 
under consideration (one-loop and two derivatives), lead to decreasing the vacuum energy density for m 2 /\eB\ > 0.927 
and to increasing that density for m 2 /\eB\ < 0.927, in accordance with Ref. S. 

Analytically, we can obtain only the limiting cases as we did in spinor electrodynamics. In particular, for m 2 -C \eB\, 
the effective action takes the following asymptotic form, 



r (2+l)scal , tj\ 
L der ( B ) 



e 2 (3 % B) 



(167r|eB|) 3 / 2 



£i (2 fc -2V2)(5-2fc)c(fc-^) 

fc=o ' L ^ ' 



2 k 



V2 



(l-2fc)C (k + 



r k 



2eB 



(94) 



while for m 2 ^> \eB\, the expansion reads 



(2+l)sco/ 



(5) 



e^B) 

' 327r 3/2 m 3 



E 



(2 2fc+3 _ 1)jB2fc+4 + (2 



2fc+l 



(2fc+l)! 



x r 2jfc 



leSl 



2/. 



l)^2feH 



(95) 



Now, let us consider the case of electric field background. Without loosing the generality, we assume that the field 
is directed along the first axis of space. We obtain 



iexp(— i7r/4) 2 



(l&n\eE\f/ 2 



I Vsinhw/ ' 



(96) 



where now uo — r\eE\ and Y = cothw. Substituting this expression into (|8l]), we come to the integral representation 
for the derivative part of the effective Lagrangian, 
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c 



(2+l)scaZ 



(E) 



ex P H7r/4)e 2 (d\\E) f duj 
(167r|e£7|) 3 / 2 



: exp 



m 
eE 



d 3 
duj 3 



d 



V sinh t 



(97) 



And we easily find the imaginary part of the expression, 



XmC 



(2+l)scal 
der 



(E) 



(d\\E) 



(-1) 



n+1 



15 + 18- 



\eE 



2 9^3| e£ ;|3/2 

1 1 n—1 



4ttV 



,5/2 



exp 



eE 



\eE\- 



TO 2 7T 3 7l 3 



eE 



m 



\eE 



which determines the correction to the corresponding result for case of constant electric field, 



lmC {2+1 ^ scal (E) 



,E\ 3/2 



o7T 

\eE\ 3 / 2 
8tt 2 



~ (_l)n+l 

„3/2 eX P 

n=l 



I) 



Li 



!3/2 



exp 



■nm 

~W\ 

„2 \ 1 



eE 



(99) 



A simple numerical analysis of the derivative correction in Eq. ( |98| ) shows that the sum in the right hand side, being 
positive for large values of the mass (or small values of the electric field), changes its sign at m 2 ps 0.721 |e£7| . Therefore, 
unlike the case of spinor QED, the time derivative of the field increases (while the gradient in space decreases) the 
probability of particle-antiparticle pair creation only for m 2 > 0.721 |e£7| . 

As in spinor QED, the two-derivative correction to the process of the pair production in scalar QED is convergent 
even in the limit of the vanishing mass. This observation, of course, is not enough to prove that the derivative 
expansion is well defined to all orders in the massless theory. 



IX. SCALAR QED IN 3+1 DIMENSIONS 



The derivative expansion for the electromagnetic field of the form (|7l]) was presented in our previous paper || (we 
just rewrite it in different form), 



(x\U bos (T)\x 



(3+1) 



1 



Itt) 2 t $ 2 sinhw " 

1 <9 M $<9^$ uj ( d 3 d 
(8tt) 2 t <P 2 [Jl! 3 + ~cLjJ 



3ujY 3 - 3Y 2 - 2ojY + l) 



J Vsinhcj/ 



(100) 



with uj = re<I>, Y — cothw. 

Now, let us consider the two most interesting particular cases as before. As in the case of the fermion theory 
presen ted i n Sec. VI, in the case of scalar QED, the derivative part of the effective Lagrangian is easily obtained by 



using (|100D with $ = iB 



r (3+i)scai , m e (d t By 

L der V B ) = 



2(87r) 2 |eB| 



dw 

— exp 



eB 



( d 3 d \ / UJ 
V duj 3 duj J V sinh u> 



(101) 



And again, as is easy to check, the situation with (101) resembles that in (2 + l)-dimcnsional scalar QED: inhomo- 
geneities of the external magnetic field lead to decreasing the vacuum energy density for large values of the ratio 

m 2 /\eB\ (m 2 /\eB\ > 0.41) and to in crea sing for small values (m 2 /\eB\ < 0.41V. 

In addition to the representation (101), we find the following one (see Eq. ( D12 ) in Appendix |d|) 



C 



(3+l)sca/ 




1 m 
^2 + 2|^ 



(102) 
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In the limit m 2 <C \eB\, this expression allows the following asymptotic expansion, 



C 



(3-\-l)scal 



(B) 



e 2 (^i?) 2 
2(87r) 2 |eB| 



18C'(-2)-ln2 



i ii 



2m 2 
3|eB| + 3|eB| 3 
2 x k 



2\eB\ 2 ^ k + 2 
1 1 fc=0 



OO 7 | 1 



1 1 fc=0 



m 



2\eB\ 

9 

m 
'2\eB\ 



(103) 



In the limit m 2 ^> le-Bl, on the other hand, we obtain 



A3 + i) SC al (m e 2 (d t B) 2 ^ (2 2fe + 3 - l)B 2k+i + (2 2 ^ 
der 1 ' ~ (87r) 2 m 2 



1)B 



2fe+2 



fe=0 



2fc+ 1 



le-Bl 



:2/r 



(104) 



In the case of the electric field directed along the first axis, we obtain 



^.(3+l)scaZ 

der 



(E) = 



ie 2 (d\\E) 2 
2(87r) 2 |e£| 



— exp 



m 

-i- -UJ 



\eE\ J V 



( dui 3 doj ) (si 



/ vsinhw 



(105) 



Thus, the imaginary part of derivative part of the Lagrangian reads 



6 + 6- 



2 7 7r 4 |e£;| 



7r 2 n 2 3 



exp 



1 



eE 



eE\ 

4 



- 1 



(106) 



which determines the correction to the probability of particle-antiparticle creation in a constant electric field expressed 
through 

lmC^ l >° al {E) = f; tl^ll exp 



16tt 3 

(eE) 



16tt 3 



2 

■Li 2 



exp 



le-Bl 



(107) 



As in (2 + l)-dimensional case, we observe that the sum in the right hand side of Eq. (106) is positive only for the 
large enough valu es of the mass (to 2 > 0.388|eS|). 

The expression (106) concludes the list of our results describing the influence of slowly varying external electromag- 
netic fields on the spinor and scalar QED vacuum in two-derivative approximation. 



X. HOW TO GET HIGHER DERIVATIVE TERMS? 



Obviously, the method of the present paper can be applied for calculating the higher derivative terms (with their 
total number equal to four or higher) of the low energy effective action in QED. However, the computational work 
with increasing the total number of derivatives is getting so hard that obtaining already all the four derivative terms 
seems to be impossible without use of a computer. Just to get feeling how difficult this problem is, let us consider 
the classification of all the relevant Feynman diagrams in four derivative approximation. 

To facilitate the calculation of the perturbative expansion in number of derivatives in the problem at hand, it is 
appropriate to develop the Feynman diagram technique. Our starting point will be the system of equations ( |l4|) and 
(jjj). Then, as is seen, the derivative expansion results from all (connected as well as disconnected) vacuum diagrams 
produced by (|l4|). A somewhat disappointing feature of our Lagrangian is an infinite number of local interactions. 
Nevertheless, as will become clear in a moment, while working at any finite order of the perturbative theory, one 
requires only a finite number of those interactions. 
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We observe that there are two different types of local interactions in ([14]). The first (bosonic) type contains only 
the bosonic fields, x^(t). The corresponding vertices are shown in Figure 1. The other interactions involve both 
the boson, x^t), and the spinor fields, t/j^t). These latter produce the vertices given in Figure 2. The integers in 
the vertices denote the number of derivatives (later called the weights of vertices) of the electromagnetic field with 
respect to space-time. Some legs in the diagrams are marked by circles and bullets. The circles correspond to legs 
related to the first Lorentz index (y) of the tensor weight, F v x tltlt ... tlin , assigned to the vertex, while the bullets, on the 
other hand, mark legs which contain the derivatives with respect to the proper time. The latter act on the (bosonic) 
propagators attached to the marked legs. 

The Feynman rules for writing expressions corresponding to Feynman diagrams are more or less standard. One 
has to use the propagators given in ( }4l| ) and ( pf2| ) for connecting the bosonic (solid) and the fermion (dashed) legs, 
respectively. The combinatoric factors can be straightforwardly derived. Two simplest diagrams giving a nonzero 
contribution to two-derivative terms in the effective action are represented in Figure 3. 

Let us mention the most general rules. To start with, we classify all diagrams leading to terms with a given 
finite number (called the weight of the corresponding diagrams from now on) of derivatives in the expansion. First 
of all, we see that diagrams with weight N may contain different number of vertices. We denote by Der(N) the 
set of all diagrams of a given weight N. By marking the bosonic (Fig. 1) and the fermion (Fig. 2) vertices with 
n derivatives just by [n] and [n], respectively, we see that the set Der(N) contains a finite number of elements: 
Der(N) = {[A], [A], [A - 1] © JlJ, [A - 1] © [1], [N - 1] © \T\, [N - 1] © [1], . . .}. Each of the elements in Der(N) 
produces in its turn a (finite) number of Feynman diagrams differing from one another by all possible connections (by 
means of propagators) between all legs of the vertices. Thus, the diagrams of weight two in Figure 3, related to C w 
and C v in the general expression (|5(]), correspond to elements [2] and [2] in the set Der(2), respectively. 

Any element of Der(N) specifics the number of different vertices as well as their separate weights. If the number of 
different vertices in a diagram is given by integers {Vi, V2, ■ . . , V*,} then the overall factor in front of the corresponding 
expression is 1/(14%! ■ ■ ■ 14!). Next, let the total number of bosonic and the fermion vertices are ks and kjr (so 
that ks + kp = k), respectively. Then the total number of bosonic legs of all the vertices in such a diagram is 
2ks + A), while the number of the fermion legs is 2kp. Since, we are interested in vacuum diagrams (with all legs 
being connected) only, the diagrams of an odd weight A are not relevant for our derivative expansion. So, we put 
A = 2n. As is easy to count, the total number of all possible connections (by means of ks + n bosonic and kp 
the fermion propagators) between these vertices is (A + 2ks — 1)!!(2&f — 1)!!, where we assume that (—1)!! = 1. 
This is an upper bound for the number of different diagrams with the given vertex set corresponding to the given 
element [Vi] © ... © [V& F ] © [I4 F +i] © ... © [Vk] € Der(N). However, due to the symmetry of the vertices with respect 
to permutations of their non-marked legs as well as with respect to permutations of identical vertices, some of the 
diagrams are in fact equivalent. For example, the naive number of all relevant diagrams for the two-derivative terms 
in the expansion of the effective action is 25. On the other hand, as is seen from our general result (|50|), the actual 
number of non-equivalent terms is 11. 

Now let us say several words about the sign factors of diagrams. First, all diagrams of weight A — 2n have an 
overall factor (—«)". To get the right sign resulting from the the fermion loops, one preliminary has to assign the 
direction of the the fermion flow in the diagram by adding arrows on the the fermion (dashed) lines. Then the overall 
sign factor is obtained by multiplying sign factors for each the fermion loop of the diagram. Each of the loop factors 
is defined by the formula: (— 1) JV »+ 1 where Ao is the number of arrows running into circles of loop vertices. This 
rule takes into account the fact that the fermion propagators are antisymmetric with respect to the simultaneous 
permutation of their Lorentz indices and proper time coordinates as well as the fact that tensor weight at the fermion 
vertices feels the order of first two indices. 

Concluding this section, we would like to express a hope that the brief description of the Feynman technique given 
here would be enough for writing a code in some of the languages used for analytical computations if such a need 
appears. 



XI. CONCLUSIONS 

In conclusion, here we further develop the method of our previous paper |J and generalize it to quantum electro- 
dynamics in 2 + 1 dimensions. The distinctive feature of our approach is the use of a special matrix basis (in Lorentz 
indices) in order to deal with functions of antisymmetric tensors such as the (background) field strength tensor in 



QED. In Sec. Ill, we give the explicit representation for these matrices as well as demonstrate how they facilitate 
the calculation. It is also the use of these matrices that allowed us to obtain the derivative expansion in the fully 
covariant form. 
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Then, in this paper, we derived explicit expression of the two-derivative term in the derivative expansion of the 
effective action in QED in both fermion and scalar QED in 2 + 1 and 3+1 dimensions. In addition, we also calculated 
the leading order corrections to the probability of the particle- antiparticle creation rate produced by space-time 
gradients of the electric field background. The latter gives a non-trivial generalization of the famous Schwinger result 
in a constant electric field ||. 

Among other results, here we derived the Feynman rules for generating the perturbative expansion of the effective 
action in the number of derivatives. This means that, in principle, an arbitrary finite order of the derivative expansion 
is calculable in our approach. For obvious reasons, the complexity of calculation explodes at higher orders and, in 
the case of the four-derivative approximation, the computational work already becomes so hard that it is almost 
impossible to get a result in the closed form without using a computer. By making use of the Feynman rules, derived 
in this paper, one can write a computer code in order to calculate higher order approximations. 

At the end, let us also make a few remarks about possible tests and applications of derivative expansion obtained 
in this paper. 

As in the case of the Euler-Heisenberg action, the derivative corrections will affect, among other things, the photon- 
photon scattering amplitude. For a vanishing background field, the latter is discussed in detail in Q. Obviously, 
when the background field is non-zero the corresponding amplitude and the energy dependence of the cross section 
are going to change. As for the explicit form of the result, it will be given elsewhere. 

Besides that, it is likely that the explicit dependence of the photon-photon cross section would be of great interest 
in studies of some real systems which exist under extremely large magnetic fields. The vicinity of the neutron stars 
and the early Universe |H| are the most natural candidates of such systems. 

The formal derivative expansion might also be useful in other problems, such as the generalization of the theory of 
magnetic catalysis of chiral symmetry breaking in QED^ |32| and QED 3 |33| to the case of inhomogeneous external 
fields. 
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APPENDIX A: COEFFICIENT FUNCTIONS WHICH APPEAR IN THE DERIVATIVE EXPANSION 



Here we give the functions 1 used in (|50j) and (p3| 



C w (a, P) = r z tanh(ar) J tf (/3r), 
C v (a, /3) = aT 3 H(aT)H(f3r) 



T 3 

Cf w (a,/3,j) = — tanh(ar)tanh(/3r)i7( 7 T), 



■[H{fir) - H{ar)], 



Cf w (a,/3,*/) = — [tanh(ar) +tanh(/3r)] 



H(ar + I3t) - H(<jt) H(<jt) 



CZ w {a,M) = — tanh(ar) ( /3tH(Pt)H(tt) 



a + f3 — 7 a + f3 J ' 

H{(3t) - Hjrrr) 
r(/3 + 7 ) 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 



1 Here we corrected the typos which appeared in ||, namely we (i) omitted an extra term in the expression for C\ v that was 
mistakenly present; (ii) replaced the wrong factor H{rf5) in the last term of (7 4 by H(tj), and (iii) added the third term in 
C^ v which was originally missing. In addition, we rewrote in a slightly different form. 
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C^i^M) = T \^ H ^ [i/i<n- //.,,) 



^ 17 (a, /9, 7) = -^H(aT)H(pr)HM - 2(/g _ ) (i?(/?r) - H{ n r)) 
r 3 (3H(Pr) r H(ar) 



2(a + 7 ) v v ' y > " 2(a + 7 )(a + /3) 
t ^(ySr) r iJ( 7 r) 

" 2(a + /3)(/3- 7 ) + 2(a + 7 )03-7)' 
yy - r 3 Q/3g(ar)g(/?T) r 3 [2(a - 0) + 7 ]g( 7 r) [/3ff(/3r) - gg(ar)] 

2 la ' P ' 7j 2 (a-/?)(a-/3 + 7 ) 2(a - /?)(a - + 7) 

olt ( 2(^ + 7) 2a-/3 + 7 



2 ^ (aT) V(a 2 -/5 2 )(a 2 -7 2 ) (a - /?) 2 (a + 7) (a - (3 \ ■ ) 
^r ff(/3r / 2(7 -a) 2/? - a - 7 



2 v " y V(a 2 -/? 2 )(/5 2 -7 2 ) (a-/3) 2 (/3-7)(a-/? + 7) 



r 



2(a/3 + 7 2 ) 7 



2 w ^( a 2_ 7 2 )(/ g 2 _ 72) (a + 7)(/3 _ 7)(a _ /3 + 7 ) 

T 



2(a-/3)(a-/3 + 7 )' 



C 3 (a, /?, 7) = — — [H((3t) - H{jt)\ 



/3 2 -~f 2 w " (a-/3)(a 2 - 7 2 ) 



/3 rr , . a/3 + 7 



2 



,., ., if(/3r) - ^ff( 7 r) 

P 7 \ a — P a ~~ 7 

, , 2ra 2 g(ar) 2r/3 2 g(/3r) 2T 7 2 g( 7 T) 

4 l a 'P'7J (a 2_ /3 2)( a 2_ 7 2) + ( a 2_ / 32)( /3 2_ 7 2) + ( a 2_ 7 2)( 7 2 _ /3 2)^ 

2 V( a + Pj( Q; +P + 7) a + 7/ 



a 1 



(a + 7 )(a + ,3 + 7 ) a + /3 



H(ar)H(-fT) 



T 



2 \a + f3 + j J [a + /?) (a + 7) 

arH(ar) (2 + g±f + g±g) r / 2g( 7 r) 
2(a + /3)(a + 7) (a + /3 + 7) + 2 ^ (a + 7 )(/3 - 7) 

2iJ(/3r) 7-ff(7T) P^pV) 



(a + /?)(/? -7) (a + 7) 2 (a + /3 + 7) (a + /3) 2 (a + [3 + 7) 
1 1 \ 



(a + /3)(a + [3 + 7) (a + 7) (a + /3 + 7) J ' 
Here we used the following notation 

. xcothx — 1 

H(x) = ^ , 

or 

and the letters with bars differ from the letters without those only in a factor of the electric charg 
that in || we ignored this difference. 

As t — > 0, these coefficient functions have the following asymptotic behavior 

C w (a, p)^^-^- (5a 2 + /3 2 ) + O (r 7 ) , 
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CY w {a,M)^^ r + 0^) 1 (A15) 

C^ w (a, M) ~ - T — + ^_ (4a 2 + 4/? 2 + 7 2 - 7a/3 - a 7 - /3 7 ) + O (r 7 ) , (A16) 

CY w (a, A 7 ) ~ ^ ( 7 - 6/3) + O (r 7 ) , (A17) 

C(M,7)--^+0(r 7 ), (A18) 

Cj V (a, 0, 7) ^ ^ ~ yjjrj ( 2 « 2 + 2 ^ 2 + 2 7 2 - 51a/3 - 9a 7 + 9/3 7 ) + O (r 7 ) , (A19) 
4t 3 t 5 

C 2 yy (a, j 8, 7) * + ^ (10a 2 + 10/3 2 + 13 7 2 + 15a/? + 3a 7 - 3/3 7 ) + O (r 7 ) , (A20) 

CX V (a, 0, 7) ^ - J + ^ (2a 2 + 2/3 2 + 2 7 2 + 9a/3) + O (r 7 ) , (A21) 

Cr V (a, p, 7) ~ £f - (a 2 + /3 2 + 7 2 ) + O (r 7 ) , (A22) 

C 5 yy (a, £,7) ~ - (20a 2 + 23/3 2 + 23 7 2 - 12a/J - 12a 7 + 6/J7) + O (r 7 ) . (A23) 



APPENDIX B: EXPANSION OF THE DERIVATIVE TERMS IN POWERS OF THE PROPER TIME 



In this appendix we give the proper time expansion of the derivative terms, as in Eqs. ( p6[ ) and (p8|), up to the 



order 



In case of spinor QED, from Eq. ( p0| ) we derive the expansion 



tr(x\U{T)\x) ~ tr(x\U{T)\x) 
ie 2 r 3 ( 7 



.,2^3 
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r r i<A,(i -r X r ufi, 



180 V 2 ' P / 315 

2F uX ^F" p p (F 2 ) x » - 2F uX ^F^P{F 2 r p - Z7F vX ^F v °>»{F 2 )\ + F vp p F ap p (F 2 ) u ' 7 



1890 



2520 



F vX , ll F a ^ p 2>%rf p F Xa F VK - Ylrf K F Xp F ila + A7ri' xp F vX F' TK + 16r/ Kp F Xa F" p 



(Bl) 



Notice that despite the difference between the two sets of matrices A'fK in 2 + 1 and 3 + 1 dimensions, the expression 
in square brackets is independent of the dimension up to this order in the expansion. In calculation, we took into 
account the Bianchi identity to show that many seemingly different terms appearing in the expansion reduce to the 
same structures. In particular, the following relations are the identities that we needed 



F \ n XK F vp = -F \ n pK F uX 



F vX>llK F^{F 2 ) XK = -F vXtlXK F v \F 2 f K , 



Fux^F^p^F^F^ = -F, J \_ IL F rTK . n n tlp F X(T F VK 



F vX ^F aK!pV VK F Xa F^ p = F^ p F aK , pV ^F Xp F^ + -F vX ^F rTK>p ri lxp F X(T F VK , 

F vX ^F aK . p ^ a F^F vX = -^F vX ^F aK , p ^ p F vX F aK , 
F vX ^F aK . p rf P F^F» x = -2F vX ^F aK . p rf p F x °F vp , 
F vX ^F a ^F^ p F vp = -jF^F^n^F^F™, 



(B2) 

(B3) 

(B4) 

(B5) 

(B6) 
(B7) 
(B8) 
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F v x,»F aK X K n ap (F 2 ) x » = -F^F^F 2 )^, (B9) 

F vX .pF aK ^ v ^{F 2 ) xp = = -lF uX ^F^>>(F 2 r p , (BIO) 

F^F aK<p r]^ K (F 2 )^ = -F uX ^F^ p (F 2 )\, (Bll) 

F^F^^^PiF 2 )^ = -F uKp F^ p (F 2 )\ + l Fu ^F^P(F 2 r p , (B12) 

F^F^p^^iF 2 )^ = F„ x ^F^P(F 2 r p , (B13) 

F^.pF aK X^ Kp {F 2 ) Xa = -F^F a /(F 2 r. (B14) 



After expanding tr (x\U (t)\x) q in Eq. (Bl) in powers of r up to the terms of order r 3 and substituting the obtained 
expression in the definition of the effective action, we arrive at the following two-derivative correction of the order 
1/m 6 , 



*-l/m« " 630m 6* W ^u\,n ' 315m 6 ^ 

j ^ F^Fr (-F" X ' fl F\ — F" X ' F V\ 

+ 270m 6 ^ *l*\2* ^ J 

„ . 2F, A , M F^ p ( J F 2 ) A ^ - 2^ A ^F^^(F 2 )^ - STF^F^F 2 )^ + F upi p F ap p {F 2 y° 



2a z 
945m 

v 2 



Fux^F^p I 387]^ F Xa F VK - l2rf K F Xp F pa + A7r) pp F vX F™ + 16r] Kp F Xa F vp ) , (B15) 



630m 6 



to the one-loop effective action in spinor QED in 3 + 1 dimensions, and the correction of the order 1/m 7 , 

-(2+1) spin _ lla 2 7T F (3~f F p u\ F ^ , jf^_l F 2\^K F u\ F 

V™ 7 " 3 36m 7 ^ + A2m 7[ ' 

+ 144m 7 ^ v 2 ^ ^' 

"2^7 ^a, m F^(^ 2 )^ - 2i^F^(F 2 r p - STF^F^F 2 )^ + F v ^F ap p {F 2 y^ 

jFvx.pFa^p ( 3Sri p ' p F Xa F VK - U^F^F^" + A7^ p F vX F aK + 16rj Kp F Xa F">* ) , (B16) 



2 

an 



336m 7 

to the effective action in 2 + f dimensions, ft turns out that these latter can be further simplified. Indeed, after 
integrating by parts, the results can be expressed through the following seven Lorentz scalars, 

L x = F^F M F» x F uX J (Bf 7) 

L-2 = F^F M F vX ^F vX ^ (B18) 

L 3 = F^F Pl F v \F vp p (B19) 

U = F vX>llK F vX {F 2 y K (B20) 

L 5 = F™F uX , p F x °F a ^ (B21) 

L 6 = F vX ^(F 3 y x (B22) 

Lt = F v ^F Xp p {F 2 Y x . (B23) 

Thus, the final results in 3 + f and in 2 + f dimensions read 

16a 2 8a 2 2a 2 , a 2 



~(3+l)spm _ — ^ — t" I — r " r 

V™ 6 315m 6 1 315m 6 2 315m 6 3 945m 6 4 

lfa 2 r 26a 2 _ 4a 2 T 

5 ~ fiTc3fi 6 ^~ 7on3fi 7 ' (B24) 



945m 6 945m 6 D f89m 6J 
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~(2+l)spm _ 2a 2 7r a 2 -K a 2 -K a 2 TT 

V™ 7 - 2 lm 6 1 21m 6 2 + 84m 6 3 504m 6 4 

lla 2 7r 13a 2 7r 5a 2 7r . . 

~504m^ 5 ~ 252m^ 6 + 126m^ 7 ' (1325) 

respectively. This should be compared with the result of || (see Eq. (14) there). Notice that the photon field in 
Ref. [|| desribes on-shell quanta, and, as a result, the terms containing L\, L$, L$ and Lj do not appear (they are 
proportional to k 2 = 0). 

In a similar way, in the case of scalar QED we obtain the following expression for the expansion of Eq. (pa) , 



(x\U bos (T)\x) ~ (x\U bos (T)\x} 



A r? n 



F vX p 



30 uX ^ 



(AF^F^ + F V \F V ^) + l -^-F vX>fiK (4T V ^F» X + 2F»\F 2 r*) 
mx,^F»P p (F 2 )^ + 13F„ x ^F» X 'P{F 2 r p + 20F V ^F^{F 2 )\ + AF vp p F Gp p [F 2 ) va 



7560 

+ l ^^F vX . p F aKiP [%rf p F Xa F VK - 4rf K 'F xp F p '° - YIrf p F uX F aK + 2Arf p F Xa F vp 
leading to the 1/m 6 correction to the effective Lagrangian density, 

2 2 

A3+i)scai _ a F 01 Fa F uX F p I a (F 2 Y R F uX F -, 
V™ 6 " 126m 6 01 ^ + 210m el ' vA ' /l * 

2 

" F^^Fr (af vX ' p f^ + F" X ' F p \ 

1080m 6 *M + A^M, J 

S^a.^V^ 2 )^ + lSF^i^OF 2 )^, + 20F uXtfi F ua ^(F 2 )\ + AF^ p F ap p (F 2 )^ 



a 2 



i/™ 7 1344m 6 1 448m 6 2 1344m 6 3 1008m 6 4 

a 2 11a 2 a 2 

"2^1^ L5+ 2016m^ L6_ 1008m^ 7 ' {BS °> 
in 3 + 1 and 2 + 1 dimensions, respectively. Here we used the same seven scalars as in the case of spinor QED above. 



3780m 6 3780m 6 1890m 6 
13a 2 a 2 r a 2 



(B26) 



3780m 6 , j 
+ 2520 6 F » X ^ F ™, P \S'n' ip F Xa F VK - Arf K F xp F^ - Yl^Pp^F ^ + 2Ar l Kp F Xa F vp \ , (B27) 
in 3 + 1 dimensions, and the 1/m 7 correction, 

-(2+1) seal _ 5a 2 7T F 0y F FV \ F fj, « 2 7T 2 x 

V™ 7 336m 7 M 112m 7i > ^ 

2 

j a 77 f 01 F,i (af vX '^F \ + F uX ' F M ^ 

"201^ ( 8 ^^^ P p(^ 2 ) Am + 13^A,^ A '"(^y p + 2QF^F^ P {F 2 )\ + AF v ^F ap p {F 2 y°^ 

~ 1344m 7 F » x ^ F ™>p {s^F^F™ - Arf~ K F Xp F^ ~ YIrf p 'F vX F™ + 247 ? Kp F Aff F 1 ''M , (B28) 
in 2 + 1 dimensions. 

Up to a divergence, the derived corrections to the effective action are equivalent to 

A3+i)scai = 13a 2 _ a 2 _ a 2 , « 2 

i/™ 6 2520m 6 1 840m 6 2 2520m 6 3 1890m 6 4 

a 2 r Ha 2 r a 2 /T ^„„ N 

^-^Zr^ + TZKKZTe 1 ^ (B29) 
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APPENDIX C: COEFFICIENT FUNCTIONS WHICH APPEAR IN PURELY ELECTRIC AND PURELY 

MAGNETIC CASES 

In this appendix we list the formulas, similar to that in Eq. (|60|), which appear in the course of reduction the general 
expression for the derivative contribution to the case of a pure magnetic (electric) field background. These are 

2 

Fvx,^ J2 C W (fjJi) A U) A (i" = - 2lCW & «) E d ^ B > ( C1 ) 

3,1 i=l 

2 

F vXtttK E CV (fi>fl) { A (f) A (i* + 2A U) A ti) = Al ° V ^ & ) E W* 3 ' ( C2 ) 

3,1 i=l 

fl^F™,, E C^UiJu fJAftAftAft = 4Cf a, a) ]T &B) 2 , (C3) 

2 

F v ^F aK , p £ CY W {h,fu fk) A (j) A (i) A (fy = -2C? w (a, -a, a) £ (^B) 2 , (C4) 

j,l } k i—l 



F vX ,,F aK , p £ Cm^fufMl-) (a$A% + A V *A* 



2 

= 2 (C™ {a, a, a) + 2C\ W (a, a, -a)) E (^ S ) 2 . ( C5 ) 

2 

F v \fj,F aK ^p £ CY w (f j ,f,J k )A%A'fiAft ) = 2C^ w (a, a, a) £ (ftS) 2 , (C6) 

i— 1 

p p \^r ,l/y ('f f, f, 1 ^4"^ 4 KCT 4 W -I- 4"^ 4 Kp 4 Aa 4- 94 yA 4 Kp 4 /io ' N \ 
^A/«,p2j i KhiJUJk) \ A {j) A (l) A {k) + A U) (i) ( fc ) + M (I) (*)J 

2 
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2 

= - (4C 2 yy (a,a,a) - C 2 vv (-a, a, a)) £ (9,B) 2 , (C8) 
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2 
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i=l 

2 
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Fy\,^F aK<p C 5 (fj,fi,fk)A(j) \ A^A^+A^A^ 
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2 

= (C 5 vy (a, -a, a) + 2C 5 yy (-a, a, a)) ]T (ftB) 2 , (Cll) 
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where a = s/IT . 

The latter expressions contain the coefficient functions from Appendix ^ calculated for a particular values of their 
arguments. The convenient representation for them reads 
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C^ v (a, —a, a 
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^vv i - - - 
G 5 (—a, a, a 

where, by definition, ui = ear and H 
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T 3 
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= — * Uco 4 H 3 + 7uj 2 H 2 - 2uj 2 H + 3H-1), 

T 3 
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APPENDIX D: SPECIAL FUNCTION REPRESENTATION FOR THE INTEGRALS WHICH APPEAR 
IN THE PURELY ELECTRIC AND PURELY MAGNETIC CASES 



In the main text, we saw that the calculation of the effective action for spinor QED in an external magnetic field 
reduces to evaluating the following integral (with yu=l/2in2 + l dimensions, and \x = in 3 + 1 dimensions) 



I^ spin \u;n) = J dujw^e-™-^ (wcothw) = J 



;/..-',,"' ' ~T~TT fcjCOth^-l-y 



oo 



I dui I ujcothu; — 1 -j — g [uj r 





d 3 



dco 3 

w\ /(3- M )(2-m)(1-m) 



uj 3 



+ 3(7(2 - M)(1 - M) +3a 2 (l- M )+^ 



(Dl) 
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where we integrated by parts (to avoid divergences as ui — > we subtracted the first two terms of the hyperbolic 
cotangent asymptotes). For large enough values of the parameter /j,, one can apply the following table integrals 



J ducj^e-™ cothuj = r(/i) 2 1 "' 1 C (ji, 1 + | 
o 

oo 

J duu^e-™ = cr-T(/i). 



(D2) 
(D3) 



Thus, the integral in Eq. (Dl) yields 

& pm \a:ji) = 2-^lV + l) 



6cr- 



1 — M 
A 6 



C M,1 + 7T 



2-/1/ CT 

- 12a -{ U-M + r 

(j, \ 2 



i k M-2.1+- 

H \ 2 



(D4) 



As one can easily check, the original integral in Eq. (Dl) is well defined for /z > — 1. Therefore, the last expression 
should allow a well defined analytical continuation to the whole that r ang e of values of \i. Notice that this should be 
true even despite the fact that the intermediate integrals, as in Eqs. (D2) and (|D3|), may not be well defined for all 
values fi > — 1. In particular, by an analytical continuation, we obtain the results for the values of y, which are of 
interest, 
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where the prime denotes the derivative of zeta function with respect to its first argument. In derivation of the second 
expression we used the following identities pTJ] 
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(D7) 



In the case of scalar QED, we come to the integral (again, with fi = 1/2 in 2 + 1 dimensions, and /i = in 3 + 1 
dimensions) 
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where we integrated by parts as in the spinor case. In addition to the table integral in (|D3j), we need also the following 
one, 
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Thus, we obtain 
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And, finally, by analytical continuation, we obtain the results for two values of /i that are of interest, 
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FIG. 1. Diagrammatic notations for the boson interaction vertices. The curly brackets denote symmetrization of the type: 

^v{X,fii . . .fj, n } = F v \^ l _,^ n + Fvy,i,\...fj, n + ■ ■ ■ + iV/z w ,fii...A- 
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FIG. 2. Diagrammatic notations for the fermion-boson interaction vertices. 
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FIG. 3. Two simplest examples of diagrams related to the two-derivative terms C w and C v in our general expression for 
spinor QED. 
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